Alternative Interpretation of Sharply Rising EO Strengths in Transitional Regions 
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It is shown that strong Of — > Of EO transitions provide a clear signature of phase transitional 
behavior in finite nuclei. Calculations using the IBA show that these transition strengths exhibit a 
dramatic and robust increase in spherical-deformed shape transition regions, that this rise matches 
well the existing data, that the predictions of these EO transitions remain large in deformed nuclei, 
and that these properties are intrinsic to the way that collecitvity and deformation develop through 
the phase transitional region in the model, arising from the specific d-boson coherence in the wave 
functions, and that they do not necessarily require the explicit mixing of normal and intruder 
configurations from different IBA spaces. 



Phase transitions are a fundamental feature of many 
physical systems and have recently been of consider- 
able interest, including extensive studies in atomic nuclei 
[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17] and in 
other mesoscopic systems, e. g. [19]. One very active area 
of study has been the study of shape changes at low en- 
ergy in nuclei [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18] 
which have been described using catastrophe [6] and Lan- 
dau theory [16, 17, 18]. Thus far, such studies have fo- 
cused on data and model comparisons for energies, E2 
transition matrix elements, and quadrupolc moments. 
However, there has been little study of EO matrix ele- 
ments in these shape transitional regions despite the fact 
that, since such shape changes are inherently linked to 
changes in nuclear shapes and radii, the EO operator and 
its transition matrix elements should provide a funda- 
mental measure of how these phase transitions proceed 
[20]. 

To study this problem, we exploit the IBA model [21], 
which provides an economic and convenient approach to 
study both phase transitional behavior, in which a simple 
two term Hamiltonian of Ising type describes transitional 
regions in terms of variations of a single control parame- 
ter, and EO transitions, for which the <i-boson content is 
explicitly related to the deformation through the intrinsic 
state formalism [21, 22, 23]. 

There have, of course, been some studies of EO transi- 
tions in transitional nuclei, most notably in the context 
of the IBA in the early work of Scholten et al. [24]. 
Their calculations for Sm isotopes provided anecdotal 
(i.e., parameter-specific) evidence for an increase in EO 
strength in deformed nuclei. Large values are also indi- 
cated in analytic expressions for ,o 2 (E0;0^ — * Of) values 
in the 0(6) and SU(3) limits [20, 21, 25, 26]. Estep et 
al [27] used calculations from ref. [28] in a shape coex- 

* ' ~ ■ P 1 

in the Mo isotopes (see below). 

However, it is the purpose of this Letter to approach 
the question of EO transitions in transitional nuclei in a 
much more general way, focusing on generic properties of 
p 2 (E0) values. We will use a simple but general Hamil- 
tonian to span the full symmetry triangle [30] of the IBA 
and will display complete contours of these monopole 



istence formalism [29] to predict p 2 (E0;0^ — > Of) values 



transitions that reveal robust, parameter- free character- 
istics of the model. As might be expected from our earlier 
comments, the most interesting behavior occurs precisely 
in shape transition regions, namely one finds a very sharp 
increase in p 2 (Ei0;0f — ► 0^ ), which then remains large 
for well-deformed nuclei. We will show that, contrary to 
common opinion, this characteristic behavior of 0^ — * Of 
E0 transition strengths does not require an explicit mix- 
ing of coexisting spherical and deformed intruder config- 
urations. Rather, it arises from a mixing of components 
with different d-boson content which is a natural ingre- 
dient in the IBA when U(5) symmetry is broken. By 
analyzing the calculated /o(E0) matrix elements in terms 
of contributions with different iid values, we will show the 
key role of the <i-boson coherence in the wave functions 
and that, while large n^ values are a necessary condition 
for large p 2 values, they are definitely not a sufficient 
condition. Finally, while surprisingly little data exists on 
Of — > Of EO transitions, much that does exist happens 
to be in shape transitional regions and we will see that 
the robust IBA predictions agree with these data. 

We start with a simple IBM-1 Hamiltonian [14] that 
includes spherical-and deformation-driving terms whose 
competition determines the resulting structure, 
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where Q = etd + dh + X (dU)^ with \ G [V7/2, 0]. 
For C — one obtains the U(5) limit while ( = 1 and \ 
= -V7/2 gives SU(3), and C = 1 and x = gives 0(6). 
In general, there is a spherical-deformed first order phase 
transition as a function of £ (except for x = where it 
is second order). The transition is most abrupt for x — 
-V7/2 and occurs at £ = 0.5 for large N, and at ( ~ 0.54 
for typical boson numbers (N ~ 10). The E0 transition 
operator is [20, 24] 



p(E0) = a(sU)° + P(dU)° = aN + p'(Sdf 



(2) 



The first term vanishes for transitions and the connec- 
tion to d-boson content is obvious. The results of our 
calculations, spanning the entire parameter space for N 
= 4, 10, and 16 are shown in Fig. 1A. 
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The essential result is immediately obvious, namely, 
that p 2 (E0;0f — + Of) rises dramatically just in the shape 
transition region, and remains large in deformed nuclei. 
This qualitative result is independent of boson number 
(i.e., N, Z), and of \- That is, there is no trajectory from 
spherical to deformed that avoids this increase. (Only the 
detailed trajectory of p 2 depends on how £ and \ vary.) 
This is a robust, parameter- free prediction of the model, 
inherent to its structure. 

These large EO transitions in the IBA raise a side is- 
sue that would be worth further exploration. The bosons 
in the IBA correspond to correlated pairs of nucleons in 
the valence space. Yet, microscopically, EO transitions 
are forbidden in a single harmonic oscillator shell [20]. 
However, realistic shell model descriptions effectively en- 
tail mixing of several oscillator shells, which is reflected 
in the use of effective charges in calculations within re- 
stricted spaces. The IBA should incorporate such effects. 

The sharp drop in p 2 for x — > at lower right (going 
toward the 0(6) limit) in the plots of Fig. 1A occurs be- 
cause of a mixing and crossing of the Of and Of states. 
This is illustrated for N =10 in Fig. IB (left). Compari- 
son with Fig. 1A shows that the Of and Of EO strengths 
interchange, and large Of — ► Of transitions emerge and 
persist into the 0(6) limit where they are the allowed 
transition from the a = (N - 2) 0+ state [20]. If the Of 
Of and Of — > Of values arc added, the contour plot re- 
mains nearly constant after the phase transition region 
[Fig. IB (right)]. Other than this case, the only strong 
ground state E0 transition is Of — > Of although strong 
transitions between pairs of excited + states abound. 

It is useful to decompose the E0 strengths in terms 
of individual components in the wave functions. This 
is done for N = 10 and x — -V? /2 in Fig. 2 which 
shows, for three £ values (one before the transition, one 
near the critical point, and £ = 1 for a well deformed 
rotor), the contributions to p(E0;0f — > Of ) from each iid 
value. These are calculated from (^(n^a^n^nd where 
<2i,2(nd) are the amplitudes in the Of and Of states with 
bosons. In U(5), the Of and 0^ states have n^ = 
and n<2 = 2, respectively and hence, by orthogonality, 
p(E0;0f — ► Of) = 0. With increasing U(5) symmetry 
breaking by £ — > 1, n^ is no longer a good quantum 
number. In fact, since £ <0^t >1 |n ( i|0f > 2 = <n 2 , > - 
<iid > 2 , the total EO strength is related to the spreading 
(fluctuations) of n^ in the ground state. As higher n^ 
components grow [14, 21, 31] so do their contributions 
to p 2 . Such rf-boson mixing is inherently related to the 
onset of quadrupolc deformation [21, 22, 23, 32]. 

Before the phase transition the p 2 values are domi- 
nated by coherent n^ = 2, 3 and 4 components. After 
the phase transition subtle positive and negative cancel- 
lations appear. Higher n^ components are essential to 
the final sum over S a\ (n^) cx2 (n<j) n^. While finite 
ci-boson amplitudes are clearly a necessary condition for 
both deformation and p values, large p values are not 
merely a trivial consequence of large <n,j > values. The 
many small p(E0;0+ -» 0+) values prove this. This is 



illustrated in the last panel of Fig. 2 which clearly shows 
the cancellations that give small p values for weak E0 
transitions. Rather, it is the specific d-boson coherence 
in the wave functions that controls the resultant p values. 

While the focus here is on Of — > Of transitions, we 
briefly comment on the behavior for higher spin. Cal- 
culations like those in Fig. 1, but for 2+ — > 2f and 
4+ — > 4f E0 transitions, show similar behavior if the p 2 
strengths are summed over all initial states. However, 
there is more fragmentation. More than one initial state 
has E0 strength to the same yrast state for a given region 
of the triangle, and different initial states dominate the 
E0 decay in different regions. Empirically, ref. [20] lists 
a number of strong E0 transitions to the first 2 + state 
and there seems to be enhanced fragmentation as well. 

The robust predictions of p 2 (E0;0f — ► Of ) demand ex- 
perimental testing. E0 Of — > Of transitions are known 
[20] in both the A = 100 and 150 transition regions. Fig- 
ure 3 compares these data with schematic IBA calcula- 
tions, using x — -V7/2 and N = 10. The data are plotted 
at £ values where the calculations reproduce the exper- 
imental R 4 / 2 values. Despite the restriction to a fixed 
X, a constant boson number N, and that £ was chosen 
simply by fitting two yrast energies, these calculations 
clearly reproduce the sharp rise in p 2 (E0) values. 

These results raise an important question relating 
to phase transitional behavior. Microscopically, the 
Fedcrman-Pittcl mechanism [33], which invokes strong 
p-n interactions [34] , leads to single particle energy shifts 
(via the monopole component [35] ) and to the descent of 
a coexisting deformed configuration in otherwise spheri- 
cal nuclei. An equilibrium deformation ensues when this 
configuration becomes the ground state. In the IBA, 
this coexistence can be explicitly included by the Duval- 
Barrett formalism [29] in which a pair of nucleons (pro- 
tons in this case) is excited across a shell or subshell gap 
to form a space with N 7rde/ = N Trsph + 2 (counting the 
extra pairs of holes and particles as additional bosons), 
thus, H = Rn b + Hat b+2 + H mix . Typical Duval-Barrett 
calculations involve many parameters - two or more for 
each term in H. The calculations of ref. [28] used 13 
parameters but reproduce the experimental p 2 (E0;0f — ► 
Of) values (see Fig. 3) in Mo rather well. 

The interesting point, however, is that, while the large 
p 2 values in these calculations have been ascribed [27] 
to the mixing, that is, to a non-vanishing H TO j x (since, 
without H m i x , E0 transitions between states of H^ =1 
and HAr 7r=3 are forbidden), it is evident from Figs. 1 and 
3 that large values of p 2 (E0;0f — > Of ) also occur in the 
IBA without the need to introduce such mixing. 

How can these seemingly conflicting results be recon- 
ciled? Figure 5 of ref. [28] shows the probabilities of 
Nn=i and components in the ground state wave 

functions (and, by orthogonality, the approximate admix- 
tures for the Of states). There is, in fact, little mixing 
(<10%) for N = 54 (spherical Mo nuclei) and even less 
(<5%) for N = 60, 62 (the first and second deformed Mo 
isotopes). Only for N = 56, 58 is there substantial mix- 
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ing. Thus, these Duval-Barrett calculations effectively go 
over into the simple (single space) IBA results before and 
after the transition region. It is therefore not the large 
p 2 value for N = 60 that requires mixing. It is rather the 
moderate p 2 values for the pre-deformed transitional Mo 
isotopes with N = 56, 58. This interpretation is validated 
by other observables. In 96 Mo 5 4, 98 Mo56, the experimen- 
tal values of the ratios B(E2;0j -» 2+)/B(E2;2+ -> 0+) 
and B(E2;2+ -> 2+)/B(E2;2+ -> 0+) exceed any pre- 
dictions of standard models, including the vibrator and 
rotor. The reason is that the Of state primarily consists 
of N bosons while the 2f , Of and 2 J states belong pri- 
marily to the N + 2 space [28] . Hence, the denominators 
are hindered. It requires the Duval-Barrett formalism 
with parameterized H m ix to account for these data. 

The key point here is that large p 2 (E0;0^ — > Of) val- 
ues in transitional nuclei can arise in two ways, either 
from mixing of coexisting spherical and intruder config- 
urations (as shown in ref. [26]) originating in different 
spaces (see ref. [36]), or, alternately, from the simpler 
IBA-1 itself. Thus, contrary to many statements in the 
literature, strong spherical-intruder state mixing is not 
required for large p 2 (E0) values, nor are large experi- 
mental p 2 (E0) values in transitional nuclei necessarily a 
signature of such mixing effects. One must analyze each 
region to determine whether to explicitly introduce shape 
mixing or whether the simple, few parameter, IBA alone 
suffices. In the Mo region ref. [28] shows that mixing of 
shape coexisting states is essential for the pre-deformed 
nuclei. However, in the first deformed nuclei in both the 
mass 100 ( 98 Sr, 100 Zr, 102 Mo) and 150 ( 152 Sm, 154 Gd) re- 
gions some of the largest known p 2 (E0) values are easily 
accounted for without such mixing by the IBA-1. 

We commented above that E0 transitions vanish in a 
single oscillator shell. It is therefore of interest to study 
how a valence space model such as the IBA, in which the 
bosons are considered to be formed from nucleons in the 
first open shell beyond an inert doubly magic core, can 
produce large E0 strengths. Of course, E0 transitions can 
arise by coupling to the giant monopole resonance, but 
this would seem to be outside the IBA space. Rather, the 
E0 transitions in the IBA may reflect the fact that realis- 
tic major shells in the independent particle model include 
an intruder orbit from the next higher shell, and that 
additional intruder orbits, from both lower and higher 
shells, appear in the Nilsson scheme with increasing de- 
formation, that is, as the phase transition proceeds. Of 
course, as a phenomenological model, one cannot relate 
the IBA directly to such a picture without detailed mi- 
croscopic analysis, but it may be that the importance of 
intruder orbits is reflected in the effective parameter, /?', 
in the E0 operator, which, in the calculations presented 
in Fig. 3, was fixed at 6 • 10~ 3 /ei?Q. Remarkably, the 



use of a simple one-body operator with constant coeffi- 
cients is sufficient for reproducing the trends of the data 
in transition regions. Given the empirical success shown 
here for the simple IBA interpretation of E0 transitions, 
microscopic studies are strongly encouraged. 

Lastly, one upshot of this study concerns well-deformed 
nuclei. The only p 2 (E0;0^ — > Of) values known in 
the deformed rare earth nuclei are very small values 
(p 2 ~ 2- 10~ 3 ) in 166 Er and 172 Yb, in contrast to the IBA 
predictions. However, the empirical Of states may not 
correspond to the Of states of the IBA, but could have 
two-quasi-particle character. Interestingly, in the neigh- 
boring nucleus 170 Yb, there is a rather strong E0 transi- 
tion (p 2 = 27(5) • 10~ 3 ) from the Of state to the ground 
state. It is also interesting that there are a number of 
large p 2 (E0;2f — > 2f) values known for deformed nuclei 
[20]. Moreover, in recent IBA calculations [37], anoma- 
lous kinks in the parameter systematics are avoided if the 
empirical O3 state is associated with the Of IBA state 
near A = 170. Clearly, it is important to measure Of — ► 
Of E0 transitions in a number of deformed nuclei to sec 
if the total E0 strength predicted in the IBA is recovered. 

To summarize, experimentally, p 2 (E0;0^ — > Of) val- 
ues rise dramatically in shape/phase spherical-deformed 
transition regions. We have presented here an alternative 
view in which this rise, and large E0 transitions in de- 
formed nuclei, arise not from mixing of coexisting spher- 
ical and deformed configurations, although such a mech- 
anism may contribute as well in specific instances (e.g., 
98 ' 100 Mo), but rather from /^-deformation and its varia- 
tion in the transition region. This result is directly con- 
nected to the physics of phase transitional regions since 
calculations within a single space reproduce the charac- 
teristic increase in E0 transition strengths. That is, using 
the IBA-1 model, we showed that contrary to common 
opinion, the rise in p 2 (E0) values is predicted even by this 
simple, single space model, that it agrees with the data, 
is parameter-free and intrinsic to the model, does not re- 
quire the mixing of different IBA spaces, and develops 
due to the specific ci-boson coherence in the wave func- 
tions. In the IBA-1 the E0 strengths are directly related 
to the fluctuations (spreading) in values, and there- 
fore to the /3-deformation. Finally, we have proposed a 
direct test of these ideas through the measurement of E0 
transitions to the ground state in well deformed nuclei. 
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FIG. 1: A) Contour plots of p 2 (E0;0^ -> 0+) throughout 
the IBA parameter space for N = 4, 10, 16. The range of 
X values implicit in the U(5) limit is explicitly shown along 
the left axis; B) Contour plots for N = 10, similar to the top 
panel, but for p 2 (E0;0^ — > 0^) on the left and for the sum 
p 2 (E0;0+ -> 0+) + p 2 (E0;0^ -> 0+) on the right. 




FIG. 2: Decomposition of the p(E0) amplitudes as a function 
of n,j for the 0^ —> Of transition for £ = 0.35, 0.54 (near the 
critical point), and ( = 1 and for 0g — > 0* for £ = 1 at lower 
right, for N = 10, x = -V7/2. Open bars are amplitudes for 
a given rid and solid bars are the running sum from n^ = 
up to the given n d value. 
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FIG. 3: Empirical p 2 (E0;0^ -> 0+) values (from ref. [20]) for 
nuclei in the A = 100 and 150 transition regions and schematic 
IBA-1 calculations. (As such schematic calculations cannot 
give R 4 / 2 <2, nuclei such as 98 Mo are not considered.) The 
solid curve is the IB A prediction for N = 10 with \ ~ ~V7 /2 
and P' (eq. 2) = 6- 10~ 3 /ei?g. The inset shows how R 4/2 itself 
behaves with (: note the similarity to the p (E0) trajectory. 
Data points are labelled with mass number A. 



